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(54) Method and apparatus for interference cancellation in a rake receiver 
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Description 
BACKGROUND 

[0001 ] Wireless communications is expanding at a phenomenal rate, as more radio spectrum becomes available for 
commercial use and as cellular phones become commonplace. For example, in the United States, wireless phone 
service is offered both in the cellular (800 MHz) and PCS (1900 MHz) bands. 

[0002] In addition, there is currently an evolution from analog communications to digital communications. Speech is 
represented by a series of bits, which are modulated and transmitted from a base station to a phone. The phone 
demodulates the received waveform to recover the bits, which are then converted back into speech. There is also a 
growing demand for data services, such as e-mail and Internet access, which require digital communications. 
[0003] There are many types of digital communications systems. Traditionally, frequency-division-multiple-access 
(FDM A) is used to divide the spectrum up into a plurality of radio channels corresponding to different carrier frequencies. 
These carriers may be further divided into time slots, referred to as time-division-multiple-access (TDM A), as is done 
in the D-AMPS, PDC, and GSM digital cellular systems. Alternatively, if the radio channel is wide enough, multiple 
users can use the same channel using spread spectrum techniques and code-division-multiple-access (CDMA). 
[0004] Direct-sequence (DS) spread-spectrum modulation is commonly used in CDMA systems, in which each in- 
formation symbol is represented by a number of "chips." Representing one symbol by many chips gives rise to "spread- 
ing," as the latter typically requires more bandwidth to transmit. The sequence of chips is referred to as the spreading 
code or signature sequence. At the receiver, the received signal is despread using a despreadingcode, which is typically 
the conjugate of the spreading code. IS-95 and J-STD-008 are examples of DS CDMA standards. 
[0005] With DS CDMA systems, coherent Rake reception is commonly used. The received signal is despread by 
correlating to the chip sequence, and the despread value is weighted by the conjugate of a channel coefficient estimate, 
removing the phase rotation of the channel and weighting the amplitude to indicate a soft or confidence value. When 
multipath propagation is present, the amplitude can vary dramatically. Multipath propagation can also lead to time 
dispersion, which causes multiple, resolvable echoes of the signal to be received. Correlators are aligned with the 
different echoes. Once the despread values have been weighted, they are summed. This weighting and summing 
operation is commonly referred to as Rake combining. 

[0006] A typical digital communications system 1 00 is shown in Figure 1 . Digital symbols are provided to transmitter 
101 , which maps the symbols into a representation appropriate for the transmission medium or channel (e.g. radio 
channel) and couples the signal to the transmission medium via antenna 102. The transmitter signal passes through 
channel 103 and is received at antenna 104. The received signal is passed to receiver 105. The receiver 105 includes 
a radio processor 106, a baseband signal processor 110, and a post processing unit 112. 

[0007] The radio processor tunes to the desired band and desired carrier frequency, then amplifies, mixes, and filters 
the signal down to baseband. At some point, the signal is sampled and quantized, ultimately providing a sequence of 
baseband received samples. Since the original radio signal has in-phase (I) and quadrature (Q) components, the 
baseband samples typically and I and Q components, giving rise to complex, baseband samples. 
[0008] The baseband processor 110 is used to detect the digital symbols that were transmitted. It may produce soft 
information as well, which gives information regarding the likelihood of the detected symbol values. 
[0009] The post processing unit 1 1 2 performs functions that depend highly on the particular communications appli- 
cation. For example, it may use the soft detected values to perform forward error correction decoding or error detection 
decoding, it may convert digital symbols into speech using a speech decoder. 

[0010] Coherent detection requires estimation of how the symbols were modified by the transmitter, channel, and/ 
or radio processor. As discussed previously, the transmission medium introduces phase and amplitude changes in 
signal, as a result of multipath propagation. The signal may also have become dispersed, giving rise to signal echoes, 
each echo having a phase and amplitude associated with it ( represented by a complex channel coefficient. Each echo 
also has a delay associated with it. Coherent demodulation requires estimation of these delays and coefficients. Typ- 
ically, the channel is modeled as discrete rays, with channel coefficients assigned to the different delays. 
[0011] A conventional baseband processor, 200, is illustrated in Figure 2. This is the standard baseband processor 
50 in a typical, coherent Rake receiver. The baseband signal is provided to a bank of correlators 202, which correlate 
different delays of the received signal to the despreading code, producing correlations, also referred to as despread 
values. The delays are provided by channel delay estimator 204, which uses known methods to estimate the delays, 
such as finding delays which give large despread values. The despread values corresponding to different delays are 
combined in combiner 206 using a weighted sum. The weights are the conjugates of channel coefficient estimates 
provided by channel coefficient estimator 208. For example, correlations to a pilot signal can be used to obtain channel 
coefficients. 

[001 2] Consider a simple example, in which the received chip-spaced baseband samples during one symbol period 
are represented by r(k). These samples are modeled as: 
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r(k) =- bc 0 s(k) + bc^s(k-1 + w(k) (1) 

where b is the symbol sent, c 0 and c 1: are the channel coefficients, the delays are 0 and 1 chip period, s(k) is the chip 
sequence used to spread the symbol and w(k) is a sequence of impairment (noise + interference) samples. 
[0013] The bank of correlators produces two despread values, denoted x 0 and x 1( corresponding to the two rays. 
These can be expressed as: 



*o = T Z>W<*) (2) 



where superscript "*" denotes complex conjugation and L is the despreading factor. Division of L is shown for illustrative 
purposes, while in practice it is well known how to extend results to the case when the division is omitted. 



1 L ' 1 

*i - T I> 1) (3) 

25 [001 4] The combiner combines the two despread values using estimates of the channel coefficients, denoted c 0 and 
a,, to produce a detection statistic that corresponds to an information symbol. This can be expressed as 
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A* A* 

z=c 0 x 0 +c 1 x 1 (4) 

The symbol value that is closest to z gives the detected value b. For BPSK modulation, bis either +1 or -1 , so that the 
detected value is given by the sign of z. 

[0015] Channel coefficients can be estimated separately using standard approaches. For example, with least mean 
square (LMS) estimation of c 0 , one would form the time varying estimate c 0 (n), where n is an index denoting symbol 
period, using 

c 0 ( n+1 )= c 0 ( n)+n S*( n)(x 0 (n)-c 0 ( n) b( n)) (5) 

A 

where ll is the LMS step size. Also b is the detected symbol value. 

[0016] It can be shown that the conventional, coherent Rake receiver is optimal when the impairment samples are 
uncorrelated. However, forcellularcommunication systems, the impairment includes interference from one's own base 
station as well as interference from other base stations. This interference is typically noise-like at the transmitter. How- 
ever, at the receiver, the interference has passed through dispersive channels, which introduce correlation. Thus, for 
cellular systems, the impairment samples are correlated and the conventional Rake receiver is no longer optimal, see, 
for example, Bottomley, "Optimizing the Rake receiver for the CDMA downlink, " Proc. 43 rd IEEE Veh. Technol. Conf. 
(VTC '93), Secaucus, NJ : May 18-20, 1993. 

[0017] Approaches which solve this problem are given in U.S. Patent No. 5,572,552 to Dentetal. Consider combining 
weight formation. First II R filtering approaches are given, in which detection statistics are formed using a weighted 
combination of despread values and a weighted combination of other detection statistics. Second, a FIR approach is 
given. Both MR and FIR approaches rely on estimating the channel responses from each base station to the receiver 
as well as noise and interference power levels. This requires multiple estimation processes that increase complexity. 
Finally, a purely adaptive scheme is given, in which the combining weights are tracked directly using decision feedback. 
However, such approaches take time to converge and don't necessarily track variations well. Thus, there is a need for 
improved combining weight computation. 

[0018] Next, consider delay estimation or correlator placement. In the aforementioned patent, an SNR criterion is 
used for tap placement that depends on channel response, noise power, and interference power estimates. Again, 
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many quantities must be estimated, increasing complexity. Thus : there is a need for a lower complexity approach to 
correlator placement. 

[0019] US 5,677, 930 discloses a spread spectrum channel estimation method having an improved method of ob- 
taining channel estimates. 

SUMMARY OF THE INVENTION 

[0020] The aforementioned problems are solved by the present invention by employing an adaptive whitening op- 
eration between the Rake fingers, so that interference that has been colored by a dispersive channel can be cancelled. 
The approach estimates fewer quantities than past approaches, lumping noise and interference into one impairment 
process, whose correlation across tap locations is estimated. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0021] The features and advantages of the invention will be understood by reading the description in conjunction 
with the drawings, in which: 

Fig. 1 is a general digital communications system; 

Fig. 2 is a baseband processor according to the prior art; 

Fig. 3 is a baseband processor according to the present invention; 

Fig. 4 is a process for correlator placement estimation according to the present invention; 

Fig. 5 is a metric computation process for use in correlator placement estimation according to the present invention; 
and 

Fig. 6 is an inhibit-controlled sliding correlator according to the present invention. 
DETAILED DESCRIPTION 

[0022] For wireless communications, the transmitter emits electromagnetic waveforms from an antenna, the medium 
is the radio propagation environment, and the receiver employs one or more antennas to recover the transmitted signal. 
While the present invention is described in the context of radio communications, it Is not limited to such systems. It is 
also applicable to wireline communications and magnetic storage systems. In such applications the radio processor 
is generalized to a device that extracts data from the transmission or storage medium. 

[0023] In the present invention , the receiver exploits the fact that, in general, the impairment (noise plus interference) 
on the different correlators is correlated. This correlation is estimated and used in the combining process. As a result, 
redundant interference components are removed. Implicitly, a whitening operation occurs, which cancels redundant 
interference components. 

[0024] A baseband processor according to the present invention, 300, is illustrated in Figure 3. Throughout, like 
reference numerals are used to refer to like elements. The baseband signal is provided to a bank of correlators 202, 
which correlates different delays of the received signal to the spreading code. The delays are provided by correlator 
placement estimator 304. The despread values corresponding to different delays are combined in modified combiner 
306. The modified combiner 306 uses channel coefficient estimates from channel coefficient estimator 208 and im- 
pairment correlation estimates from impairment correlation estimator 310 to combine the despread values to form a 
detected symbol output. 

[0025] The impairment correlation estimator 31 0 estimates correlation between the impairment on the different cor- 
relator outputs. In the past, impairment correlation among different antenna signals has been used with MLSE reception 
to cancel interference, as disclosed in U.S. Patent No. 5,680,419 to Bottomley. 

[0026] With the present invention, correlation among different despread values from the same antenna signal is 
exploited in a Rake receiver. 

[0027] The approach for estimating the impairment correlation can be similar to the approaches given in the Bottomley 
patent, except that despread values from the same antenna are used instead of received samples from different an- 
tennas. For example, error signals for the different correlators can be used. Forthe two-ray example, these error signals 
would be: 
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= x 1 -bc 1 



(7) 



where b is a detected symbol value. If there are pilot symbols, known values can be used. If there is a pilot channel, 
there is an effective symbol value, usually + 1 . 

[0028] Collecting the error signals into a vector e(n), where n denotes symbol period, then an impairment correlation 
matrix estimate R(n)can be updated using: 
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R(n)=XR(n-1)+e(n)e n (n) 



(8) 



where superscript "H" denotes Hermitian transpose. Note that R(n) is Hermitian, i.e. R H =R, so that only the diagonal 
and one of the off diagonal triangles (upper or lower) need be estimated and updated. For the two-ray example, the R 
matrix has the form: 
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R = 



Poo Poi 

POI Pi! 



(9) 



25 



where p 00 and p u are real numbers (imaginary part is zero). Note that the inverse of this matrix is used in combining. 
The inverse can be updated directly using the well known matrix inversion lemma. The term "impairment correlation" 
is used to refer not only to correlations, but to any related quantities, such as an inverse correlation matrix. 
[0029] The modified combiner 306 would then combine the despread values using both the channel co-efficients 
and the impairment correlations to produce a detection statistic. The combining operation can be expressed as: 



30 



:S H R- 1 X 



(10) 



35 



where the channel coefficients have also been collected into a vector. The detection statistic can be further processed 
to produce a detected symbol value. It can also be used as a soft value for further processing. 
[0030] The combining can be implemented in a number of ways. The despread values can be combined by the 
impairment correlation matrix first, then combined by the channel coefficients. Alternatively, the channel coefficients 
and the impairment correlation can be pre-combined to form weights: 



£-1 a 
w=R c 



(11) 
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Then, the detection statistic can be expressed as: 



z'=w"x 



(12) 



Alternatively, combining with the impairment correlation matrix estimate and the channel coefficient estimates can be 
done together. For the two-ray example, this can be expressed as: 



i^j t(Pii c o-Poi c i)*o + (Poo<VPoi c o)*i] 



(13) 



where 
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IR| -PooPn-|Poi| (14) 
[0031] Also, it is possible to factorize the inverse impairment correlation matrix estimate using square-root factori- 
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zation, giving two factors: 
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Thus, the combining operation can be expressed as: 
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where 



y=Qx 
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and 
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g=Qc 



(18) 
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Multiplying the despread vector by Q whitens the impairment, but changes the overall channel response. As a result, 
the channel coefficients have to be modified as well. 

[0032] Based on square-root Kalman filtering, it is possible to estimate and track the square-root matrix O. which is 
another form of impairment correlation. Then, the despread values used for detection and the despread values used 
for channel coefficient estimation can be whitened first. Standard channel coefficient estimation applied to the whitened 
despread values will yield g. Alternatively, one can track the channel, then apply the square root to the channel coef- 
ficient estimates, as shown above. 

[0033] The operation of the correlator placement unit 304 can employ any conventional delay estimation approach. 
For example, unit 304 can employ the approaches disclosed in pending U .S. Patent Application Serial No. 09/005,580, 
entitled "Multiple Delay Estimation for Direct Sequence Spread Spectrum Systems" filed on January 12, 1998 
[0034] However, the correlator placement unit 304 can alternatively employ a different approach, based on modifying 
the approaches in the aforementioned pending application, e.g. : modifying the metric to include an estimate of the 
impairment correlation matrix. 

[0035] For example, one of the approaches is illustrated in Figure 4. The process starts in start block 402. Then 
despread values corresponding to different delays are generated and stored in step 404. Then, at step 406, hypothe- 
sized tap locations or delays are made. In step 408, a metric is then calculated corresponding to this set of delays. In 
step 410, the metric is compared to previous ones. If the metric is better, it is stored as the new best metric and the 
corresponding delay combination is also stored. Then, in step 412, it is determined whether all the delay combinations 
have been exhausted. If not, the next combination is considered in step 406. Otherwise, the delay estimates are taken 
to be the ones corresponding to the best metric and the process ends in step 418. 

[0036] The key distinction is how the metric is calculated. The metric calculation step 408 is further detailed in Figure 
5. For the delay combination, a set of chan nel coefficients are estimated using standard approaches in step 502. These 
channel coefficients should be the "composite" channel coefficients, corresponding to the transmit, medium, and re- 
ceive filter responses. Side information in the form of transmit and/or receive filter response knowledge can be used 
to improve estimation, as discussed in the aforementioned pending U.S. Patent Application to Sourour et al. Then, in 
step 504, the impairment correlation is estimated using the approaches described previously. Finally, in step 506, the 
metric is computed using both the channel coefficient estimates and the impairment correlation estimates. The preferred 
metric can be expressed as: 



which corresponds to an SNR figure of merit. 

[0037] Similarly, impairment correlation among different antenna signals has been used with synchronization, as 
disclosed in pending application U.S. Patent Application Serial No. 08/773,560 to Bottomley and Chennakeshu, filed 
December 27, 1996. 

[0038] With the present invention, correlation among different despread values form the same antenna signal is 
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exploited to determine correlator or "finger" placement in a Rake receiver. 

[0039] It may be advantageous to perform conventional delay estimation first, then consider delays within a certain 
proximity of the conventional delay estimates. It may also be desirable to keep the M strongest ray delay estimates, 
then only consider alternatives for the remaining P delay estimates. 

5 [0040] Referring back to Figure 3, the bank of correlators 202 can be realized in a number of ways. It can be a group 
of integrate-and-dump correlators. It can also be realized using a single sliding correlator. In this case, delays associated 
with the bank of correlators correspond to selecting which of the outputs of the sliding correlator are kept for further 
processing. A third approach is to use a selectively inhibited sliding correlator as illustrated in Figure 6. 
[0041] Data samples are provided to delay line 602, which includes delay elements 604a-604c. For this example, it 

10 is assumed that the sampling rate is two samples per chip and that the despreading code section has length three. It 
will be apparent to one skilled in the art that this aspect of the present invention may be extended to any sampling rate 
and any despreading length. Also note that delay 604a may be omitted, depending on how the input samples are 
generated. 

[0042] The delayed samples are provided to processing engine 606, which includes remove chip units 608a-608c 
15 and adder 610. The delayed samples are provided to remove chip units 608a-608c, where the despreading chip values 

are removed from the samples, producing modified samples. For example, the received sample is multiplied by the 

conjugate of the despreading chip value. When chip values are +1 or -1 , then chip removal is simply negating or not 

negating the received sample. The modified samples are added together in adder 610, to produce despread values. 

[0043] What distinguishes the operation of the processing engine 606 from a conventional sliding correlator is that 
20 the operation of the elements in the processing engine 606 can be selectively inhibited, saving power when despread 

values are not needed. Thus, remove chip units 608a-608c and adder 610 have a control input, which determines 

whether an operation will be performed or not. 

[0044] The processing engine 606 is controlled by inhibit unit 612, which produces the control signal based on the 
delays to be used. In essence, the inhibit control unit 612 instructs the processing engine 606 to produce despread 
25 values only for those delays of interest. For all other delays, the inhibit control unit 61 2 instructs the processing engine 
606 not to process the delayed samples. 

[0045] The delay line 602 can be efficiently implemented as a circular buffer. This avoids the power consuming need 
to repeatedly shift data samples. 

[0046] This approach overcomes the problem that a bank of 4 integrate-and-dump correlators can only examine up 
30 to 4 delays. It also overcomes the problem that a sliding correlator has, which is large power consumption. Theinhibit- 
controlled sliding correlator can be used with conventional Rake combining as well as for initial acquisition and delay 
tracking. 

[0047] Referring back to Figure 3, other forms of correlation combining can be used, based on known antenna array 
processing approaches. For example, the impairment correlation estimate can be replaced by a despread correlation 
55 estimate, in which the error signal e is replaced with x when forming the R matrix. This approach will also cancel 
interference, though the "soft" detection statistic will not work as well in subsequent processing as the preferred em- 
bodiment already given. 

[0048] The present invention can use any type of channel coefficient tracking algorithm. For example, the LMS, 
KLMS ( see , e.g., Jama/ etal., "Adaptive MLSE performance on the D-AMPS 1900 channel," IEEE Trans. Veh. Technot., 

40 vol. 46, pp. 634-641 , August 1997), RLS, and Kalman tracking algorithms are appropriate. While chip-spaced rays 
were used as examples, the rays can have arbitrary spacing , includingf ractional spacing . Channel coefficient estimation 
can also be done using interpolation between pilot symbol sections. Similarly, the present invention can use a number 
of approaches to impairment correlation estimation. The impairment correlation can be either tracked or interpolated 
between pilot symbol sections. When receiver quantities are tracked; per survivor processing (PSP) can be used to 

45 improve performance, by keeping channel coefficient estimates and impairment correlation estimates per hypothesized 
symbol values. 

[0049] A modified approach is possible, in which the impairment correlation estimate used to combine the despread 
values is a combination of the impairment correlation estimate and a fixed value. This provides a way of gracefully 
switching between conventional approaches (the fixed value is the identity matrix) and the present invention. It can 

50 also be used to switch between an adaptive estimate and a known structure. For example, if the interference is non- 
dispersive, then the impairment is colored only by the receive filter. Thus, the fixed matrix could be a matrix of pulse 
shape autocorrelation values, possibly scaled by an estimate of the noise power. The "fixed matrix" could also be 
adaptive, using a pulse shape autocorrelation matrix scaled by an adaptive noise power estimate. 
[0050] The present invention can be used in a multi-pass approach. Despread values corresponding to a data frame 

55 can be stored. In the post-processing phase, forward-error-correction and forward-error-detection decoding can be 
used to correct or detect errors. Then, re-encoding can be used to provide reference symbols for a second-pass, for 
better parameter estimation. Multi-pass demodulation is described in U.S. Patent No. 5,673,291 to Dent which is in- 
corporated herein by reference. 
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[0051] The present invention can also be used in conjunction with multiple receive antennas. In pending U.S. Patent 
Application Serial No. 08/992,174 to Jonas Karlsson and Sara Mazur, entitled "Code division multiple access mobile 
station interference suppression", multiple antenna despread values corresponding to a particular delay are combined 
using an impairment estimate across antennas. With the present invention, the despread values from all antennas 
would be collected together into one set of despread values, which would be combined according to the present in- 
vention. Thus, impairment correlation across antennas and across delays would be estimated and used in combining. 
[0052] A hybrid approach is possible, in which groups of despread values are combined using the present invention, 
and those groups are then simply added together to form the detection statistic. Unlike the aforementioned application 
to Karlsson et al., the groups do not necessarily have to correspond to the same delay, but different antennas. 
[0053] The invention has been described in the context of a single, modulated traffic channel. However, the invention 
is also applicable to systems with pilot symbols or with a pilot channel, such as the IS-95 downlink. With pilot symbols, 
the symbol values are known, so that channel coefficient estimation and impairment estimation can use known, instead 
of detected, symbol values. With a pilot channel, the pilot channel can be viewed as a continuous sequence of known 
symbol values (usually all +1). Thus, these known symbol values could be used. 

Claims 

1. A method for despreading a received signal comprising the steps of: 

storing a plurality of data samples; 

selectively processing the stored data samples by combining them with chip values; and 
summing the result to produce despread values. 

2. The method of claim 1 wherein said step of selectively processing further comprises the step of: 

selectively combining said stored data samples with chip values during each sample period. 

3. The method of claim 1 , wherein said step of selectively processing further comprises the step of: 

selectively combining said stored data samples only for delays of interest. 

4. A sliding correlator comprising: 

a delay line having a plurality of delay elements for delaying a stream of data samples input thereto; 

a plurality of remove chip units, each associated with an output of one of said plurality of delay elements, for 

removing chips from said delayed data samples; 

an adder for receiving outputs from said plurality of remove chip units and adding said outputs together; and 
an inhibit unit for selectively enabling at least one of said plurality of remove chip units and said adder. ' 

5. The sliding correlator of claim 4, wherein said delay line is a circular buffer. 

6. The sliding correlator of claim 4, wherein said inhibit unit operates to enable said at least one of said plurality of 
remove chip units and said adder only for delays of interest. 

7. The sliding correlator of claim 4, wherein said inhibit unit selectively enables and disables said plurality of remove 
chip units and said adder together. 
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Fig. 4 
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Fig. 5 
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